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Abstract. The present work explores the theoretical effects of rotation in calculating the period ratios of double-mode radial 
pulsating stars with special emphasis on high-amplitude d Scuti stars (HADS). Diagrams showing these period ratios vs. periods 
of the fundamental radial mode have been employed as a good tracer of non-solar metallicities and are known as Petersen 
diagrams (PD). In this paper we consider the effect of moderate rotation on both evolutionary models and oscillation frequencies 
and we show that such effects cannot be completely neglected as it has been done until now. In particular it is found that even for 
low-to-moderate rotational velocities (15-50kms"'), differences in period ratios of some hundredths can be found. The main 
consequence is therefore the confusion scenario generated when trying to fit the metallicity of a given star using this diagram 
without a previous knowledge of its rotational velocity. 

Key words. Stars: variables: S Set - Stars: rotation - Stars: variables: RR Lyr - Stars: oscillations - Stars: fundamental 
parameters - Stars: variables: Cepheids 



1. Introduction 

Radial pulsators, and more particularly double mode pulsators, 
have been extensively studied using the well known Petersen 
diagrams (from now on called PD). Such diagrams show the 
ratio between the fundamental radial mode and the first over- 
tone as a function of the fundamental mode. Typically, the stars 
concerned are double-mode Cepheids, RR Lyrae and high- 
amplitude 6 Scuti stars (HADS). Recently, the analysis of data 
from large-scale proj ects like OGLE ( Optical Gravi tational 
Lensing Experiment ISzvmanskil l2005t JUdalski et alJ fT9 97-). 
NSVS (Northern Sky Variability Survey Wozniak et al. 2004), 
AS AS (A ll Sky Automated S urvey Poima nskL2002i..2003.) or 
MACHO llAlcock et alJl2000h has permitted deeper studies of 
the observational properties of such double-mode pulsat ors. 
These period ratios were firstly studied by .Petersen> d 19731 



Il978h . and they have been used for decades as metallicity in- 
dicators, as a function of stellar mass and age, to test mass- 
luminosity and/or radius-luminosity relations. In addition, they 
also have been used, for instance, to determine the distance 
modulus to the SMC (Kovacs 2000). 

While non-radial pulsators like /3 Cephei or low-amplitude 
6 Scuti stars (LADS) generally show moderate to fast rotational 
velocities (vsin/), the double-mode radial pulsators can be con- 
sidered as slow-to-moderate rotating stars. This fact has lead 
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to neglect systematically the effect of rotation on theoretical 
period-ratios, in particular F-IO ratios in PD. However, we re- 
mind the reader that such objects are generally faint enough 
as to make difficult the observation of their rotational veloci- 
ties (vsin/). Figure [2 shows the location of all known double- 
mode HADS (from Poretti et al. 2005| in the HR diagram, 
compared with the location of all LADS with measured vsin/. 
As can be seen, the narrow band (in effective temperature ) 
occupied by double-mode HADS overlaps the region where 
LADS are located. In luminosity, HADS occupy a broader 
range. Unfortunately, the number of known HADS constitutes 
a very poor sample, specially when compared with LADS. 
This means that we cannot discard, a priori, the possibility of 
HADS in a wider range of effective temperatures and present 
larger vsinj measurements. In addition to this, although most of 
HADS present vsin/ < 20kms this represent a lower limit 
of their rotational velocities. When varying the angle of incli- 
nation of the star /, velocities up to 50 km s ' could be reached. 

From the theoretical side, only a few works examined par- 
tially the effect of rotation on period ratios and ma i nly fo cused 
on very rapidly rotators. In lPerez H ernandez et al.l(ll995l) . sec- 
ond order effects of rotation on oscillation periods were taken 
into account in order to discriminate between radial and non- 
radial modes. The authors expressed the relative period change 
as 5nJ,'°'Vn„ - Z„(ni-ot)/n„)"^, where D,, is the unperturbed 
period of the mode with radial order n, and n(rot) being the ro- 
tation period. The coefficient Z„ depends on the radial order 
of the mode and on the structure of the star. Such values were 
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Fig. 1. Absolute magnitudes My as a function of the effective 
temperature of all known S Scuti stars (empty symbols), com- 
pared with those of all double-mode HADS known up to date 
(filled symbols). Differe nt symbol type s represent the observed 
vsin/ ranges (taken from lRodrfguez et al-2000.) . The cross rep- 
resents typical errors on absolute magnitudes and effective tem- 
perature for 6 Scuti stars. 

estimated interpolating from comput ations deri ved for a poly- 
tropic model with index 3. Recentlv. iPamvatnvkh (.2003 ) stud- 
ied the behaviour of period ratios of radial modes when near 
degeneracy effects due to rotation are included for a typical 
1 .8 Mq 6 Scuti stellar model. He showed that very large and 

non-regular pe rturbations to such ra tios are expected to occur. 
I 1 J 1 

More recently, Suarez et al. r2005a) proposed a limit of valid- 
ity of the perturbation theory (up to second order) in terms of 
rotational velocity for rotating models. Such limit is given by 
the behaviour of period ratios when near degeneracy is consid- 
ered, which clearly complicates the naive interpretation of the 
PD. 

In the present work, we aim at analysing the consequences 
of neglecting the effect of rotation on radial period ratios 
even for low rotational velocities. To do so, up-to-date tech- 
niques taking properly into account the rotation in the mod- 
elling (in equilibrium mo dels and in oscillation fre quencies) are 
used. Si milarly as done in lPamvatnvkhl ( I2OO3I) and lSuarez et alJ 
ll2005ah near degeneracy is also considered. 

The paper is structured as follows: A general description of 
the modelling, focusing on how rotation is taken into account is 
given in Sect. |2l Fundamental-to-first harmonic ratios in pres- 
ence of rotation, IIi/o {Q), are introduced in Sect. |3] and a dis- 
cussion of their impact on PD analysis is proposed in Sect. |4] 
Finally, the conclusions are given in Sect.|5] 

2. The modelling 

The evolutionary code CESAM ("MorefigQ?^ is used, which 
is particularly adapted for our purposes. The numerical preci- 
sion and the mesh grid (around 2000 mesh points given in the 
basis of B-splines) of equilibrium models have been adapted 
according to the oscillation computation requirements. 



Fig. 2. Typical PD (Ilo in d) containing different tracks of 
I.8M0 evolutionary models computed with different initial 
metal content [Fe/H], from solar composition (bottom, solid 
line) to -1.00 (top, triple dot-dashed line). Crosses represent 
non-rotating models. The shaded area corresponds to typical 
values found for Pop. I stars. The two filled star symbols rep- 
resent the observed Ili/o of the double-mode high-amplitude 
6 Scuti stars RV Ari and SX Phe, as an illustration of (Pop. I) 
and (Pop. II) stars respectively. (For clarity, colours are used in 
the on-line version of the paper). 

Following iKippenhahn & WeigertI ( fl990l) . a first order ef- 
fect of rotation is taken into account in equilibrium equations. 
In particular, the spherical symmetric contribution of the cen- 
trifugal acceleration is included by means of an effective grav- 
ity guS - g - ^c{r), where g represents the local gravity com- 
ponent, r the radial distance and Jicir) - | Q, the centrifugal 
acceleration of matter elements at a distance r from the centre 
of the star. This spherically symmetric contribution of the ro- 
tation does not change the shape of the hydrostatic equilibrium 
equation. Although, the non-spheric components of the cen- 
trifugal acceleration are not considered, they are included as 
a perturbation in the oscillation computation. The total angu- 
lar momentum of models is assumed to be globally conserved 
along the evolution of the star. Input physics has been adapted 
for intermediate mass stars. 

Theoretical oscillation spectra are computed from the equi- 
librium models described in the previous section. For this 
purpose the oscillation code Filou (Tran Minh & Leon 1995; 
ISuarezl2002l) is used. This code, based on a perturbative anal- 
ysis, provides adiabatic oscillations corrected for the effects of 
rotation up to second order (centrifugal and Coriolis forces). 

Furthermore, for moderate-high rotational velocities, the 
effects of near d egeneracy are expected to be significant 
(iSoufi et al.ll998l) . Two or more modes, close in frequency, are 
rendered degenerate by rotation under certain conditions, cor- 
responding to selection rules. In particular these rules select 
modes w ith the same azim uthal order m and degrees I differ- 
ing by 2 dSoufi et al.lll998l) . If we consider two generic modes 
a\ = (n,{,m) and 0-2 = in' ,1' ,m') under the aforementioned 
conditions, near degeneracy occurs for Ictq-j - ctqj | < ctq, where 
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cr„| and cr„, represent the eigenfrequency associated to modes 
ai and 0-2 respectively, and ctq represents the stellar rotational 
frequency. In certain cases, such effect may be dom inant in the 
behaviour of th e radial period ratios studied here ( P amvatnvkhl 
l2003USuarez" et al. 2005a). However, due to its complexity, we 
believe such effects should be analysed separately (Suarez et 
al., in prep.) and, for the sake of clarity, they have not been 
included in the present work. 

3. The rii/o (Q.) ratios 

In general, low order radial period ratios can be considered as 
dependent of the distribution of mass (or density) throughout 
the star, and the thermodynamical properties of the stellar mat- 
ter. In Fig.|2j a classic PD for 1.8 evolutionary tracks is de- 
picted. It illustrates the well-known dependence on the metal- 
licity of I Ti/o = n i /On r a tios similar as those shown for in- 
stance in Petersen ( 1973), Petersen & Christensen-Dals gaardI 
(Il996l) and Petersen & Christensen-Dalsgaard (.1999) . As can 
be seen, the period ratios increase when decreasing the stellar 
initial metal content. This property is commonly used to dis- 
criminate, in the context of radial pulsators. Pop. I from Pop. II 
stars. The shaded region indicates the typical Ili/o values found 
for Pop. I stars, in the range of IIi/o = [0.772, 0.776]. As a ref- 
erence, the observed period ratios of two stars: RV Ari (Pop. I) 
and SX Phe (Pop. II) are also depicted (values obtained from 
IPoretti et al. 2005. and references therein). For a given chem- 
ical composition the period ratio IIi/o is determined by mass 
and the radius (or Do which is scaled by the mean density p, 
and then 

Hi/o = ni/o(M,/?,Z). 

In addition, period ratios can be easily written in terms of 
observed quantities through the so-called pulsation constant 
Q„^( - n„ f ^Jp/pQ■ In the context of radial modes, this con- 
stant can be expressed as 



allowing us to assign to each observed radial period n„ a value 
Qnfi, which can be compared with theoretical predictions. 

However, it is well known that rotation modifies the struc- 
ture of stars and thereby the cavity where modes propagate. It is 
thus plausible to consider the period ratios to be Q-dependent. 
Moreover, the angle of inclination of the star is generally un- 
known (only rotational projected velocities (vsin/) are provided 
by observations), and thereby it follows 

ni/o(n) = ni/o(M,«,z,Q(0). 

From a theoretical po int of view, following ISoufi etalJ lll99}j): 
LSuarez et al.l ll2005bl) . the adiabatic oscillation eigenfrequen- 
cies cJn/^m can be expressed in terms of a perturbative theory 



_ (0) , (2) 



for radial modes (f - 0, m - 0), where y represents the 

unperturbed frequency and <wjf(,(, the second order correction 
term. Both terms implicitly include the symmetric component 
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Fig. 3. Fundamental radial order mode n (top panel) and the 
first overtone uji ci (bottom panel) normalised frequencies as a 
function of the logarithm of the fundamental radial order period 
(in d). (For clarity, colours are used in the on-line version of the 
paper). 



of the centrifugal force (mainly the departure from sphericity of 
the star) given by the equilibrium model. First order correcting 
terms w''^ „ (Coriolis force effect) are proportional to the az- 
imuthal order m, and therefore they are zero for radial modes. 
For shortness, from now, only the subscript n is kept. Figure|3] 
illustrates the effect of rotation on normalised' frequencies cor- 
responding to the fundamental mode wo.n (top panel) and the 
first overtone a>i,n (bottom panel) for a given metallicity. As 
can be seen, the first overtone frequencies are visibly more af- 
fected by rotation than the fundamental ones. Without enter- 
ing in details (a theoretical work analysing the behaviour of 
ITi/o (fi) is currently in preparation) this suggests the possibil- 
ity that W] n is proportional to the distribution of the density in- 
side the star (affected by rotation) rather than simply the mean 
density. 



' The normalisation constant is (GM/i?^)''^ which is proportional 
to the model mean density p. 
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Fig. 4. Theoretical PD including rotation effects (RPD). The IIi/o period ratios have been computed in the manner described 
in Sect. 121 for a set of evolutionary 1.8, M© tracks obtained for different metallicities. Tracks for three initial rotational veloc- 
ities are considered: 25 and 50kms"' (from top to bottom). Left panels show only rotating models. Right panels show the 
comparison between rotating and non-rotating tracks (classic PD). For convenience, the following symbols are used: crosses, 
representing non-rotating models; asterisks, representing models evolved with Qi - 25 km s ' ; diamonds and those evolved with 
Qi - 50kms As in Fig|2] filled star symbols represent the observed values for SX Phe and RV Ari. (For clarity, colours are 
used in the on-line version of the paper). 



4. The ITi/o (fl) period ratios and the metallicity 
determinations 

In order to analyse the impact of considering IIi/o (f2) rather 
than rii/o period ratios on metallicity determinations, these 
quantities, obtained from models taking into account rotation 
Hi/o (Q), are compared with those Ili/o. To do so, several evo- 
lutionary tracks have been computed as described in Sect.|2]for 
six different metallicities: [Fe/H] = 0, -0.1, -0.2, -0.35, -0.50 
and -l.OOdex; and three different initial rotational velocities 
Qi = 25,50 and 100 kms '. The case of rotational velocities 
between 50 and 100 km s"' is purely illustrative, since there are 
no HADS known in that range. The discussion is thus mainly 
focused on models with rotational velocities up to 50, km s ' . 

During the evolution, the rotational velocity of models de- 
crease up to 0.75 f2i at TAMS, due to the global conservation of 
the total angular momentum (see Sect.|2|for more details). The 



mass of models is fixed to 1 .8, M© which typically corresponds 
to a (5 Scuti star. 

Adiabatic oscillations are then computed from these rotat- 
ing models, obtaining thus the corresponding ni/o(0) period 
ratios. In Fig. |3 such rotational PD (hereafter RPD) are dis- 
played, from top to bottom, for tracks computed from Qj - 25 
to 100 kms"' respectively. It can be noticed that ni/()(Q) pe- 
riod ratios increase for increasing rotational velocities (left pan- 
els). Such effect is similar to decreasing the metallicity in clas- 
sic PD. As explained in the previous section, this can be under- 
stood in terms of variations of the density distribution in stellar 
model interiors due to rotation effects (Suarez et al., work in 
preparation). In addition, the shift to larger period ratios is de- 
pendent on the metallicity. In particular, the higher is the rota- 
tional velocity the closer are the models of different metallicity 
in RPD (left panels). This means that the effect of rotation on 
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Fig. 5. RPD (Ho in d) illustrating the effect of considering dif- 
ferent masses. Symbols have the same meaning of those from 
Fig. 13 except for empty and filled circles which correspond 1.9 
and 2.00 M0 non-rotating models and solar metallicity respec- 
tively. (For clarity, colours are used in the on-line version of the 
paper). 

period ratios is systematically larger for increasing metalUcity 
values (up to the solar value, in the present study). 

In Fig.l^Ili /o (il) (right panels) are displayed together with 
classic ones Ili/o. A first quantitative comparison is performed, 
in terms differences of period ratios 



(5n,/o(Q, [Fe/H]) 



Hi/o (O) - n 



1/0 



[Fe/H] 



for a given metallicity. For the lowest initial rotational velocity 
considered, Qj - 25kms"', OTi/o reach up 2 - 3 lO""', and for 
Qi - 50kms"', differences increase up to 6 - 8 10"^ 

, and finally, for Q; = lOOkms they can be of the order 
of 10"^. The main effect on PD (right panels) is to shift and 
compress tracks of same metallicity (and mass) toward higher 
period ratios with respect to classic tracks. 

The previous differences in period ratios can also be anal- 
ysed in terms of metallicity. For shortness, tracks will be spec- 
ified, from now on, with the subscript corresponding to the ro- 
tational velocity considered. For instance, the track computed 
with Qi = 25kms ' and [Fe/H] - -0.1 will be called as 
[-0.1]25. As it will be shown, such analysis leads to a confus- 
ing scenario: For Q; = 25kms"' (top, right) solar metallicity 
models tracks are similar to classic ones. However, when de- 
creasing [Fe/H], rotating and non-rotating tracks are located 
quite close. In particular, rotating [-0.10]25 tracks may be 
confused with [-0.20]o ones, and [-0.20]25 may be confused 
with non-rotating [-0.35]o ones. The analysis of the middle 
right panel reveals that [0.00]5o tracks are located closely to 
[-0.20]o ones. Similarly, [-0.10]5o tracks may be confused 
with [-0.20, -0.35]o ones, and finally, [-0.20]5o tracks are 
close to [-0.50]o ones. As can be seen, the mix-up is critic for 
pop. I stars when considering 1 .8 rotating models evolved 
withQi = 25, 50 km s'. 

Up to this point, all the discussion has been based on the re- 
sults obtained only for 1 .8 Mq models. When other masses and 



metallicites are taken into account, the confusion in metallici- 
ties and rotational velocities significantly increases. In order to 
illustrate this, two solar metallicity tracks for higher mass (1.9 
and 2.00 Mq) non-rotating models are displayed in the RPD of 
Fig. |5l where they are compared with other rotating non-solar 
tracks. As can be seen, the effect of increasing the mass of the 
models goes in the same direction as increasing the rotational 
velocity, that is, it increases the period ratios. It is worth not- 
ing that such shift toward larger period ratios systematically 
occurs for any rotating track. Therefore the previous discus- 
sion based on 1.8 M© models is equivalent whatever the mass 
and/or rotational velocity considered. Such behaviour may ex- 
tend the confusing scenario to Pop. II stars. Considering 1 .8 Mq 
models, it would be necessary to consider initial rotational ve- 
locities up to 100 kms '. However, this lower limit rapidly de- 
crease when increasing the mass of the models. For instance, 
2Mo and Qi - 50 kms ' tracks may be misinterpreted with 
[-1.00]o ones. 

The construction of complete RPD for a wide range of 
metallicities, initial rotational velocities and masses becomes 
necessary for the exhaustive analysis of double-mode pulsators 
(work currently in progress), however this exceeds the scope of 
this paper. 

5. Conclusions 

The impact of taking into account the effect of rotation on 
Petersen Diagrams has been examined here, focusing on main 
sequence double-mode pulsators. Detailed seismic models 
have been computed considering rotation effects on both equi- 
librium models and on adiabatic oscillation frequencies. For 
1.8 M0 stellar models, period ratios have been calculated for 
different rotational velocities (Rotational Petersen Diagrams) 
and metallicities, and then compared with classic non-rotating 
ones (PD). 

The analysis of these RPD reveals that the difference in 
period ratios increases with the rotational velocity for a given 
metallicity. It remains around lO^^* for rotational velocities up 
to 50 kms"' and it can reach up 10"^ for rotational velocities 
close to 100 km s ' . Such difference have been found enough to 
produce a significant confusing scenario when analysing RPD 
in terms of metallicity variations. In particular, for 1.8 Mq stel- 
lar models, differences in metallicity up to (5[Fe/H] ~ 0.30 dex 
can be found when considering models evolved with initial ro- 
tational velocities of 50 kms"'. Furthermore, such confusion 
may still increase when including other stellar masses, rota- 
tional velocities and metallicities, as well as other physical 
parametrisation (work in progress). 

The results presented here should thus be taken into ac- 
count when analysing double-mode pulsators with Petersen 
Diagrams, in particular when accurate metallicity and/or mass 
determinations are required. 

A work on the detailed analysis of period ratios in pres- 
ence of near degeneracy is currently in preparation. Such work 
would provide new constraints on the modelling of HADS. 
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